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ABSTRACT 

o 

Photometric reverberation mapping employs a wide bandpass to measure the AGN continuum variations and a suitable band, usually a 
narrow band (NB), to trace the echo of an emission line in the broad line region (BLR). The narrow band catches both the emission line 
, and the underlying continuum, and one needs to extract the pure emission line light curve. We performed a test on two local AGNs, 

X^/^ ' PG0003+199 (=Mrk335) and Arkl20, observing well-sampled broad- (B, V) and narrow-band light curves with the robotic 15cm 

telescope VYSOS-6 on Cerro Armazones, Chile. In PG0003+199, Ha dominates the flux in the NB by 85%, allowing us to measure 
the time lag of Ho- against B without the need to correct for the continuum contribution. In Arkl20, H/3 contributes only 50% to the 
flux in the NB. The cross correlation of the B and NB light curves shows two distinct peaks of similar strength, one at lag zero from 
the autocorrelated continuum and one from the emission line at T cmt = 47.5 ± 3.4 days. We constructed a synthetic H/3 light curve, by 
subtracting a scaled V light curve, which traces the continuum, from the NB light curve. The cross correlation of this synthetic HB light 
curve with the B light curve shows only one major peak at T ce „, = 48.0 ± 3.3 days, while the peak from the autocorrelated continuum 
at lag zero is absent. We conclude that, as long as the emission line contributes at least 50% to the bandpass, the pure emission line 
i-C . light curve can be reconstructed from photometric monitoring data so that the time lag can be measured. For both objects the lags 

we find are consistent with spectroscopic reverberation results. While the dense sampling (median 2 days) enables us to determine 
JL T cmt with small (~10%) formal errors, we caution that gaps in the light curves may lead to much larger systematic uncertainties. We 

calculated a virial black hole mass for both objects in agreement with literature values, by combining ^blr (jcent) from photometric 
monitoring with the velocity dispersion of a single contemporaneous spectrum. We applied the flux variation gradient method to 
estimate the host galaxy contribution in the apertures used and thus the host-subtracted restframe 5 100 A luminosity L AGN during the 
time of our monitoring campaign. At the AGN/host contrast of our sources (L AGN / 'Li, os , ~ 1) the uncertainty of L AGN is ~10%. For 
both sources, our L AGN differs significantly from previous estimates, placing both sources ~50% closer to the /?blr _ L AGN relation. 
Our results suggest that the scatter in the current Rblr - L AGN relation is largely caused by uncertainties in Rblr due to undersampled 
light curves and by uncertainties in the host-subtracted AGN luminosities inferred so far. If the scatter can be reduced sufficiently, 
then two quasar samples matching in Rblr should match also in intrinsic L AG n, independent of redshift, thus offering the prospect 
to probe cosmological models. Photometric reverberation mapping opens the door to efficiently measure hundreds of BLR sizes and 
host-subtracted AGN luminosities even with small telescopes, but also routinely with upcoming large survey telescopes like the LSST 
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1 . Introduction Reverberation based BLR sizes and black hole masses have 

been determined so far for only about 45 low- and intermediate- 
luminous AGN (Wandel et al. 1999, Kaspi et al. 2000, Peterson 
Reverberation mapping (Blandford & McKee 1982) has revolu- e t al. 2004, Bentz et al. 2010). Six years of reverberation map- 
tionised our understanding of active galactic nuclei (AGN) dur- ping of six luminous high-redshift (z > 2) quasars, using the 
c3 ; ing the past decades. Spectroscopic monitoring measures the 9.2 m Hobby-Eberly-Telescope in Texas, resulted in a tentative 
mean time lag r between brightness changes of the triggering R BhR determination of only one quasar, while in the other five 
optical-ultraviolet continuum (from the compact accretion disk sources the existence of pronounced variation patterns could be 
around the black hole) and the emission line brightness echo established (Kaspi et al. 2007). 
of the BLR gas clouds further out; the mean BLR size is then 
Rblr = t-c, where c is the speed of light. The BLR size as mea- 
sured by the centroid T ce „, of the cross correlation of light curves Reverberation spectra usually require observations with at 
represents a luminosity-weighted average (Koratkar & Gaskell least a 2m-class telescope. The BLR sizes range between a 
1991, Penston 1991). Evidence exist that the BLR clouds are few and several hundred light days and the time spans become 
essentially in virialised motion around the central black hole even longer in the observer's frame due to the time dilation fac- 
(Peterson et al. 2004 and references therein). Combining Rblr tor 1 + z. To match distinct echo features, light curves over a 
with velocity dispersion 5v of the emission line, the virial black few months to several years are needed. Consequently, spectro- 
hole mass can be estimated as Mbh = / ■ ^blr ■ Sv/G, where / is scopic reverberation mapping is resource expensive and it be- 
a scaling factor in the order of unity to account for the geometry comes prohibitive at high redshift. Therefore, more efficient ob- 
of the BLR and G is the gravitational constant. serving strategies are desired. 
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Another issue with respect to the BLR size - AGN luminos- 
ity relationship is how to determine the AGN luminosity free of 
host galaxy contributions. While high spatial resolution images 
have been invoked to model and subtract the host contribution 
from the total fluxes (e.g. Bentz et al. 2009), the flux variation 
gradient method (Choloniewski 1981) provides an efficient ap- 
proach which we will inspect here as well. 

We here revisit photometric reverberation mapping. It em- 
ploys wide pass bands to trace the AGN continuum and suit- 
able narrow-bands to trace the echo of the BLR lines. The ad- 
vantage is that photometric monitoring can be obtained much 
faster and with small telescopes. Furthermore, while long-slit 
spectroscopic monitoring allows for simultaneous measurement 
of typically one or two calibration stars placed accurately on 
the spectrograph slit, one may expect that the large number of 
non-variable calibration stars on the same images and of similar 
brightness as the AGN facilitates to attain photometric monitor- 
ing with high precision. 

About 40 years ago, narrow-band AGN monitoring of emis- 
sion lines has been performed, for instance for NGC3516 and 
NGC4151 (Cherepashchuk & Lyutyi 1973). Despite the coarse 
(>10 days) time sampling, these observations discovered a lag 
between Ha and the continuum of 25 and 30 days for NGC 35 16 
and NGC 4151, respectively. However, spectroscopic monitor- 
ing with better time sampling finds a lag of 13 ±4 days 
for NGC 3516 (Peterson et al. 2004) and 6.6 + 1.1 days for 
NGC 4151 (Bentz et al. 2006). The discrepancy between the 
photometrically and spectroscopically determined lags is most 
likely due to insufficient time sampling of the early photometric 
monitoring. 

To explore the capability and achievable accuracy of pho- 
tometric reverberation mapping, we selected two local AGNs, 
the relatively clear-cut case PG0003+199 (Mrk335) at redshift 
Z = 0.0258 and the more challenging case Ark 120 at z = 0.0327. 
Spectroscopic reverberation results are available for comparison 
(Peterson et al. 1998a, 2004). 

2. Data 

We performed a monitoring campaign between August 2009 
and March 2010 using the robotic 15 cm telescope VYSOS-6 
of the Universitatssternwarte Bochum on Cerro Armazones, the 
future location of the ESO Extremely Large Telescope in Chile. 
VYSOS-6 is equipped with a 4096 x 4096 pixel CCD yielding a 
field-of-view of 2.7°, and seven broad- and narrow-band filters. 1 

We obtained light curves with a median sampling of 2 days 
in the fi-band (Johnson, 4330 + 500 A), the redshifted Ha (NB 
6721 + 30 A at z = 0.0258) and UB lines (NB 5007 + 30 A at 
z = 0.0327). For Ark 120 we also obtained a light curve in 
the V-band (Johnson, 5500 ± 500 A). In addition, we observed 
PG0003+199 in both B and V on July 25' h 2010 and June 26* 
201 1 . FigureQ]shows the effective transmission of the filters used 
here. To perform absolute photometric calibration, each night 
we observed standard stars in the fields SA092, SA095, SA111 
from Landoldt (2009). For both targets, a single contemporane- 
ous spectrum was obtained with CAFOS at the 2.2 m telescope 
on Calar Alto, Spain, with a slit width of 1'.'54. 

We reduced all data in a standard manner, using IRAF and 
custom written tools. Because the flux calibration using the stan- 
dard star fields introduces additional errors in the light curves, 
we created relative light curves (in normalised flux units) using 

1 Since August 2010, VYSOS-6 consists of two 15 cm telescopes on 
the same mount, equipped with 14 broad- and narrow-band filters. 
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Fig. 1. Effective transmission of the filters convolved with the 
quantum efficiency of the ALTA U16M CCD camera. 
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Fig. 2. Observed spectrum of PG0003+199. The B-band 
(blue shaded) traces essentially the AGN continuum, while the 
narrow -band (red shaded) is dominated by the Ha line. Note that 
the blue and red shaded areas are a schematic illustration and 
that quantitative calculations were made using the transmission 
curves shown in Figure Q] 



20-30 non-variable stars located on the same images within 30' 
around the AGN and of similar brightness as the AGN. For the 
analysis of time lags we used the mean and standard deviation of 
these relative light curves. For the photometric analysis (to ob- 
tain the AGN luminosities), we kept the shape of the mean light 
curves fixed and calibrated them by a least-squares fit to the pho- 
tometry derived from the standard star fields. Source parameters 
and photometry results are summarized in Table [T] 



3. Results 

3.1. Results on PG0003+199 

PG0003 + 199 is essentially a point-like low-luminosity quasar 
(Narrow Line Seyfert-1) without a bright extended host. The Ha 
line is strong and well covered by the narrow-band filter. This 
makes PG0003+199 a clear-cut test case. 
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Fig. 3. B-band and Ha light curves of PG0003+199. The Ha 
observations were started in October 2009 after recognizing the 
pronounced B-band variations. 



3.1.1. Spectrum 

Fig. [2] shows the spectrum of PG0003+199. The contribution of 
higher order Balmer lines (and that of the host galaxy) to the 
B-band is negligible, thus the B-band is dominated by the AGN 
continuum. Taking into account the spectral resolution, the ob- 
served Ha line dispersion <x = 1300km/s reduces to an intrinsic 
value cr = 870km/s. The narrow -band filter effectively covers 
the line between velocities -2800 km/s and +1800km/s, so that 
at least 95% of the line flux is contained in the band pass, as we 
determined after line profile deconvolution. 

The contribution of both the [N II] 6583 A and narrow-line 
Ha flux is predicted to be smaller than 40% of the [O III] 
5007 A emission (Bennert et al. 2006), hence negligible (< 
10%). The continuum underneath the emission line is small, 
contributing to only 15% of the total flux in the band pass. 
Therefore, the narrow-band light curves will be dominated by 
the Ha echo of the BLR gas clouds. Comparison with simulated 
line profiles of echo models (Welsh & Home 1991, Home et al. 
2004) ascertains that the 5% line flux outside the band pass has 
only a marginal effect (< 2%) on the BLR size determination: 
The narrow-band echo may miss only a small fraction (< 20%) 
of the innermost part of the BLR, namely that part which exhibits 
the fastest line-of-sight velocity, while the innermost gas clouds 
moving closer along the sky plane exhibit a modest line-of-sight 
velocity and are therefore contained in the narrow-band. 



3.1 .2. Light curves and BLR size 

Fig. |3]depicts the light curves of PG0003+199. The AGN con- 
tinuum as traced by the B-band increases gradually from August 
2009 to a peak around end of September, followed by a steep 
trough by 20% at MJD -55118 (= 55045 + 73). Also the Ha 
light curve has a steep ~10% trough with a delay of about 20 
days. 

We used the discrete correlation function (DCF, Edelson & 
Krolik 1988) to cross correlate the Ha and B-band light curves. 
The cross correlation shows a major peak with a lag of 20.2 days 
as defined by the centroid T cen , (Fig. @). Two smaller correla- 
tion peaks around lag 50 days and 65 days are present; they are 
caused by the B-band troughs at begin and mid of September 
2009 (Fig. [3). We here do not consider these two lags further, 
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Fig. 4. Cross correlation of the B and Ha light curves of 
PG0003 + 199, taking all data points. The dotted lines indicate the 
error range (±lcr) around the cross correlation. The red shaded 
area markes the range to calculate the centroid T cent of the major 
correlation feature, yielding a lag of 20.2 days (vertical dashed 
line). 



because spectroscopic reverberation predicts a lag of less than 
30 days (Peterson et al. 2004). 

To determine the lag uncertainty, we applied the flux ran- 
domization and random subset selection method (FR/RSS, 
Peterson et al. 1998b). From the observed light curves we cre- 
ated 2000 randomly selected subset light curves, each containing 
63% of the original data points, and randomly altering the flux 
value of each data point consistent with its (normal-distributed) 
measurement error. Then we cross correlated the 2000 pairs of 
subset light curves and computed the centroid T cent . Fig.|5]shows 
the distribution of the 2000 r cent values which yields the median 
lag r cent and the 68% confidence range. From this distribution 
we obtain a lag T cent = 20.5 days. Thus, the mean of pos- 
itive and negative uncertainty of T cen , is about 12% (2. 4/20. 5). 2 
Correcting for the time dilation (factor 1 +z at z = 0.0258) yields 



a rest frame lag of 20.0 ^'7 days. 

We did not find Ha lags of PG0003+199 in the literature, 
but we can compare our result with that of the H/3 line for 
which spectroscopic reverberation mapping found a rest-frame 
lag of 15.7 ± 3.7 days averaged over several epochs (Peterson 
et al. 2004, Bentz et al. 2009). This H/3 lag is -30% smaller 
than our Ha lag of 20.0 days. One possible explanation could 
be that the BLR size is smaller for higher excitation emission 
lines. Kaspi et al. (2000) found in the spectroscopic reverbera- 
tion data of 17 PG quasars Ru„ = s ■ Bh/j with a scaling factor 
s = 1.19, but a high uncertainty of +0.23. On the other hand, 
Kollatschny (2003a) provides clear evidence for the BLR strati- 
fication of MrkllO (s = 1.37, his Table 1). Recently, the Lick 
AGN Monitoring Program of 11 low-luminosity AGN found 
s = 1.54 + 0.4 (Bentz et al. 2010). It could be that the BLR 



2 The observed Ho- light curve has only 10 data points, three of which 
are in the trough. While a subset of exactly 63% data points would re- 
quire to take 6.3 data points, we selected 7 random data points. This 
~10% larger number of data points may bias the errors to somewhat 
smaller values. On the other hand, we find the largest errors for those 
randomly selected subset light curves, where all three data points in the 
trough are omitted. In this case the information about the existence of 
a trough is completely lost. To counterbalance this effect, we selected 7 
instead of 6.3 data points. 
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Fig. 5. Result of the FR/RSS error analysis of the cross corre- 
lation of the B and Hor light curves of PG0003+199. The his- 
togram shows the distribution of T cent for 2000 flux randomised 
and randomly selected subset light curves. The median of this 
distribution is taken as the final value of the lag r cent . The red 
shaded area marks the 68% confidence range used to calculate 
the errors of T cent . 



stratification and the scaling factor s depends on the luminosity. 
The luminosity of PG0003+199 is similar to that of Mrk 1 10. If 
the value s = 1.37 of Mrk 110 holds for PG0003+199, then the 
HB lag of PG0003+199 translates to an Hor lag of 15.7 x 1.37 = 
21.5 days, hence agrees within 8% with the restframe Ha lag of 
20.0 days from our measurement. 

3.1.3. Black hole mass 

Next we try to determine the mass of the central black hole. 
While spectroscopic reverberation data enable one to use the 
velocity information from the rms spectra, i.e. from the vari- 
able portion of the emission lines, several attempts have been 
performed to determine the mass of the central black hole us- 
ing single epoch spectra (e.g. Vestergaard 2002, Woo et al. 
2007, Denney et al. 2009a). Thereby the reported uncertain- 
ties of the velocity dispersion range between 15% and 25%. 
Combining for PG0003+199 the intrinsic line dispersion from 
our single spectrum (870km/s, with 25% uncertainty adopted) 
and the time delay (20.0 + 2.4 days) yields a virial black hole 
mass Mbh = 2.8 + 1.1 • 1O 6 M , in agreement with the value of 
2.8 + 0.8 • 10 6 M o derived via spectroscopic reverberation map- 
ping by Peterson et al. (2004). 

3.1.4. Host-subtracted AGN luminosity 

To determine the AGN luminosity free of host galaxy contri- 
butions, we applied the flux variation gradient (FVG) method 
proposed by Choloniewski (1981). In this method the B and V 
data points obtained in the same night through the same aper- 
tures are plotted in a fi-flux versus V-flux diagram. Note that 
fluxes (e.g. in units of mJy) are plotted and not magitudes. The 
important feature is that the flux variations follow a linear re- 
lation with a slope Ybv given by the host-free AGN continuum 
(e.g. Choloniewski 1981, Winkler et al. 1992, Winkler 1997). 

In the flux-flux diagram the host galaxy - including the con- 
tribution of line emission from the narrow line region (NLR) - 
lies on the AGN slope somewhere toward its fainter end. This has 



McAlary15.7" " 




V flux [ mJy ] 

Fig. 6. B versus V fluxes of PG0003+199, measured in a 7'.' 8 
aperture (black dots). Errors are within the size of the symbols. 
A linear least-squares fit to the three data points yields the range 
for the AGN slope plotted by the two steep lines. This slope 
is corroborated by the slope (dashed line) obtained with larger 
apertures 14" (Elvis) and 15'.'7 (McAlary). The range of host 
slopes plotted by the two shallow lines is taken from Sakata 
et al. (2010). The intersection area of the AGN and host slopes 
(red shaded) defines the host contribution. The blue shaded area 
marks the intersection of the AGN slope with the range of B 
fluxes during our Aug-Nov 2009 monitoring campaign. This 
area allows us to infer the range of V fluxes during the mon- 
itoring campaign. The vertical dotted line marks the host flux 
found in an aperture of 5" x 7'.' 6 via modelling of HST F550M 
images by Bentz et al. (2009). All data are corrected for galactic 
foreground extinction. 



been demonstrated for a sample of 1 1 nearby AGN (neither con- 
taining PG0003+199 nor Arkl20) by Sakata et al. (2010), using 
monitoring data in combination with Hubble Space Telescope 
and multi-band ground-based images from the MAGNUM tele- 
scope. For 8'.'3 aperture, they derived host galaxy colours in the 
range 0.8< B—V < 1.1 (their Table 6), which agree with the typ- 
ical colours of a bulge or an elliptical galaxy. This colour range 
corresponds to a host slope in the range 0.4 < T h °y < 0.53. The 
contribution of the NLR emission lines to the B and V bands is 
less than 10% of the host flux and has a similar color as the host 
galaxy (cf. Tables 6 and 8 of Sakata et al. 2010). We note that 
dust enshrouded nuclear starbursts, which often accompany the 
AGN phenomenon, have red colors (B — V ~ 1). They would 
be unresolved on HST images and could explain the small de- 
viations of the host position from the AGN slope found in three 
cases by Sakata et al. (2010). Therefore, by host we here denote 
the host galaxy including NLR AGN lines and starbursts. 

The host slopes pass through the origin. Because the host 
slopes are flatter than typical AGN slopes (Vbv ~ 1), the inter- 
section of the two slopes should occur in a well defined range. 
Winkler et al. (1992) showed for NGC3783 that this intersec- 
tion range allows one to determine the host flux contribution and 
to calculate a host-subtracted AGN luminosity at the time of the 
monitoring campaign - even without the need for high spatial 
resolution images. Next (and in Sect. 13.2.4b we test this proce- 
dure on our data. Table [2] summarizes the results of the FVG 
diagnostics, which are derived as follows: 
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During our monitoring campaign of PG0003+199 in Aug- 
Nov 2009, we obtained only B fluxes, but we have three data 
points at hand for the B versus V flux-flux diagram as shown in 
Fig. [6] All data are corrected for galactic foreground extinction 
(Tab.Q). The fact that PG0003+199 was much brighter in 1979 
(McAlary et al. 1983) than in 2009-201 1 allows us to achieve a 
linear least-squares fit to the data points with sufficient accuracy. 
The fit yields an AGN slope F BV = 1.20 + 0.04. This slope is at 
the steep end, but within the range determined by Winkler et al. 
( 1 992) for other type- 1 AGN. 

To check the slope of PG0003+199 further, we searched in 
the NASA Extragalactic Database (NED) and found two more 
B and V data points taken with comparable apertures, one with 
15'.'7 by McAlary et al. (1983) and one with 14'.'0 by Elvis 
et al. (1994). These (extinction corrected) data points are dis- 
placed to the right of the 7'.'8 slope in Fig. [6] because they con- 
tain additional host emission from the extended disk which is 
redder than the AGN emission. These two data points yield a 
slope Fbv = 1 -25 ± 0.04. The consistency of Fbv derived from 
two independent data sets corrobborates the steep AGN slope of 
PG0003+199. 

The range of host slopes (0.4 < F' B °y < 0.53) is taken from 
Sakata et al. (2010). The intersection area of the AGN and host 
slopes defines the host flux. Averaging over the intersection area 
yields a mean host flux of 2.16 + 0.29 mJy in B and 4.56 + 0.30 
mJy in V. 

During our monitoring campaign, the total B fluxes lie in 
the range between 5.69 mJy and 6.78 mJy with a mean of 6.24 + 
0.3 1 mJy. To derive the contemporaneous 5 100A AGN luminos- 
ity, which is widely used in reverberation studies, we tried to ex- 
trapolate the host-subtracted B fluxes to restframe 5 100 A fluxes, 
using a power-law spectral shape (F v oc v°) with a being con- 
strained by Fbv- But this extrapolation results in large uncertain- 
ties, because a depends sensitively on errors in Fbv- Therefore 
we applied the following procedure, inferring in a first step the 
V flux range. 

During our monitoring campaign, the V fluxes should lie in 
the intersection area of the B flux range with the AGN slope 
(blue shaded area in Fig. |6j. The fact that the data point from 
June 2011 lies in this area supports our suggestion. From this 
area we infer the V fluxes during the monitoring campaign, 
yielding a mean fV - 7.95 + 0.31 mJy. 

Then, during our monitoring campaign the host-subtracted 
AGN fluxes are in the range 4.08 + 0.43 mJy in B and 
3.39 + 0.40 mJy in V. From this range we interpolate 
the host-subtracted AGN flux of PG0003+199 at restframe 
5100A f 510 oA - 3.52 + 0.41 mJy, adopting for the interpola- 
tion that the AGN has a power law spectral energy distribution 
(F v oc v ff ). To determine the errors on F 5100A , we interpolated 
■^5iooA between fB + cr and fV + cr as well as between fB - cr 
and fV - cr. At the distance of 112.6 Mpc this yields a host- 
subtracted AGN luminosity i 5100 A = 3.13 ± 0.36 x I0 43 erg/s. 
Note that the 12% uncertainty includes the measurement errors, 
the uncertainty of the AGN and host slopes, and the AGN varia- 
tions. 

The AGN luminosity during our campaign derived with 
the FVG method is about a factor two smaller than the value 
L 5100A = 6.03 + 0.03 ■ 10 43 erg/ s derived with host galaxy mod- 
elling by Bentz et al. (2009) for the data of the spectroscopic 
monitoring campaign (Peterson et al. 1998a). To understand this 
difference, we compared the numbers given in Bentz et al. (their 
Tables 7-9) with ours. Our aperture area (7'.' 8 in diameter) is 
only 25% larger than that of Peterson's campaign (5" x7'.'6), and 



the additional flux from the outer region is negligible (Fig. 3 in 
Bentz et al. 2009). The total observed (not extinction corrected) 
fluxes ^ 5I00A are similar, being 7.83 x ergs' 1 cmr 2 A~ l 
for our campaign and 7.68 - 8.81 x 1Q' 15 ergs' 1 cm' 2 A' 1 for 
Peterson's campaign. The main difference, however, lies in the 
estimate of the host contribution. Already in the early stage of 
the HST image decomposition with GALFIT, Bentz et al.'s host 
flux (fV host ~ F55QM = 1.88 x 10~ 15 ergs' 1 cm' 2 A" 1 ) corre- 
sponds after extinction correction to 2.12m/>', which is more 
than a factor two smaller than our fVhost - 4.56mJy. 

As illustrated in Fig. [6] it is hard to bring the small host flux 
from GALFIT modelling into agreement with the available FVG 
data. Firstly, we consider whether the discrepancy may be caused 
by failures of the FVG method. We find: 1) Altering the assump- 
tions on the host slope does not bring the host contributions into 
agreement. 2) Although the AGN slope is determined by only 
a few data points, the large flux range, which allows for a good 
slope determination, and the consistency between the AGN slope 
determined with apertures of 7'.'8 and ~15" argues in favour of 
a steep AGN slope. 3) The correction for galactic foreground 
extinction may introduce errors in the AGN slope. However, the 
foreground extinction is relatively small (Table[T} and even with- 
out correcting for foreground extinction the AGN slope remains 
steep (Fbv = 1-16 ± 0.05) resulting in fVhost = 4.2 + Q.3mJy, 
hence about twice the GALFIT host value. 4) If the host flux 
from GALFIT modelling were correct, then at low luminosity 
the AGN slope (of PG0003+199) would show a strong curvature 
towards redder colours, in contradiction to the results by Sakata 
et al. (2010) where the host galaxy lies on the linear extension of 
the AGN slope. These four arguments lead us to conclude that 
the discrepancy between GALFIT and FVG host fluxes is not 
caused by errors in the FVG method alone. On the other hand, 
looking for potential error sources in the GALFIT method, it 
could be that a large fraction of the bulge and/or nuclear star- 
bursts of PG0003+199 is unresolved on the HST images and 
therefore underestimated by the GALFIT modelling. Future data 
may clarify this issue. 

3.2. Results on Ark 120 

3.2.1 . Implications from the spectrum 

Arkl20 lies at redshift z = 0.0327 so that the HB line falls into 
the NB 5007+30 A filter. Compared to the case of PG0003+199, 
the observations of Ark 120 are more challenging: Firstly, com- 
pared to Ha, the HB line is fainter contributing only to ~50% 
of the flux in the NB band pass (Fig. |7). In order to measure 
the continuum variations underneath the HB line, light curves 
using neighbouring off-line intermediate-bands would be ideal 
but such bands were not available. While the V-band covers HB 
and the Fe complex, it is dominated to 83% + 3% by the AGN 
and host continuum (calculated using the transmission curve of 
filter and CCD shown in Fig. Q]) and therefore we used the V- 
band, in addition to B-band and NB 5007 A. The second chal- 
lenge is that (after accounting for the spectral resolution) the HB 
line is broader than the narrow bandpass, so that the broadest line 
wings, carrying about 10 - 15% of the line flux, escape detection 
in the NB. Despite of these handicaps, the results are promising. 

3.2.2. Light curves and BLR size 

Fig. [8] depicts the light curves of Ark 120. The dense time sam- 
pling clearly outlines the variations. The B-band steeply in- 
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Fig. 7. Observed spectrum of Ark 120. The B-band (blue 
shaded) traces the AGN (and host) continuum and a small (< 
5%) contribution by Hy. The V-band (green shaded) traces the 
continuum and a small (< 15%) contribution by Hp. The narrow- 
band traces the continuum (black shaded) and Hp (red shaded), 
each component contributing to about 50% to the band pass. The 
red wing of Hp extending beyond [OIII] is neither seen in Ho- 
nor in Hy and could be due to Fe lines (Korista 1992). For the 
effective transmission curves of B and V see Figure Q] 
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Fig. 8. Light curves of Ark 120 between Oct. 2009 and March 
2010. The B, V and NB curves are as observed, the HB curve 
is computed by subtracting a scaled V curve from the NB curve 
(HB = NB - 0.5 V) and re-normalization (mean=l). For better 
visibility, the light curves are shifted against each other by +0.2. 
The December 2009 data are mostly lost due to strong wind pre- 
venting us from opening the dome. 



creases during ten days November 2009 by about 17% with a 
sharp peak at end of November 2009; from January to March 
2010 it declines gradually by about 17%. In November the V 
and NB light curves are similar to that of the B-band, but with 
smaller amplitude (V ~ 10%, NB ~ 7%). In contrast to the steep 
B band flux increase, the NB flux increase is stretched until 
January 2010 with an amplitude of 12-15% between begin of 



November and January. The possible influence of the missing 
data in December 2009 is discussed in Section|4] 

The NB contains about 50% continuum and 50% HB (Fig.[7Ji. 
In order to remove the continuum, we computed a synthetic HB 
light curve by subtracting a scaled V curve from the NB curve: 
HB - NB - 0.5 V (V scaling factor = 0.5). Alternating the V scal- 
ing factors between 0.4 and 0.6 yields similar results. The syn- 
thetic Hp curve still has a remaining variability pattern in Oct- 
Nov which we suggest to arise from the relicts of the continuum. 
However, the important result is that compared to November the 
~10% higher flux level in January becomes more pronounced. 
While the subtraction of a scaled V curve from the NB curve 
turns out to be a successful approach to recover a pure Hp light 
curve, we discuss potential refinements in Section|4] 

Fig. [9] shows the cross correlation results applying in a first 
step the DCF to the entire light curves. The DCF of B and NB 
exhibits two distinct major peaks, one at a lag of 47.5 days and 
one at a lag of about zero days. Because the cross correlation of 
B and V yields a lag of 1 .2*1 g days consistent with lag zero, we 
conclude that the lag zero peak in the DCF of B and NB comes 
from the continuum contribution contained in the NB. This con- 
tribution appears here quasi in autocorrelation, while the lag of 
47.5 days is due to the line echo. This is further corroborated by 
the DCF of the B and synthetic Hp curves, where the peak at lag 
47.7 days dominates, while the continuum caused contribution 
at lag zero is much reduced. This shows that the subtraction of 
a scaled V light curve from the narrow-band light curve, in fact, 
largely removes the interferring continuum emission, allowing 
us to measure the lag of the emission line. 

To determine the lag uncertainties, we applied the FR/RSS 
method, again creating 2000 randomly selected subset light 
curves, each containing 63% of the original data points, and ran- 
domly altering the flux value of each data point consistent with 
its (normal-distributed) measurement error. Then we cross cor- 
related the 2000 pairs of subset light curves and computed the 
centroid T cent . As shown in Fig. [TO] the resulting median lag is 
T ce „, = 47.5 +3'o days for the DCF of B and NB, and T cenl = 48.0 
+ Yl days for the DCF of B and synthetic Hp. For short, the lag 
uncertainty is about 7%. Correcting 48.0 days for the time di- 
lation factor (1.0327) we obtain a rest frame lag of 46.5 + 3.25 
days, consistent with the lag found by spectroscopic monitoring 
(47 ± 10 days by Peterson et al. 1998a, 2004, and in the range 
34-54 days by Doroshenko et al. 1999). 

3.2.3. Black hole mass 

From our single epoch spectrum we determined an intrinsic Hp 
line dispersion cr = 1950km/s, in agreement with the value 
by Peterson et al. (2004). To get rid of the red wing of the 
Hp line, we mirrored the blue side of the Hp line profile to 
the red, and calculated the line dispersion of this synthetic pro- 
file. The shape and the dispersion of this profile are consis- 
tent with what we derived for Ha {cr ~ 1970 km/s) and Hy 
(cr ~ 1920 km/s). Combining line dispersion (with 25% uncer- 
tainty adopted) and the time delay (46.5 + 3.25 days) yields a 
virial M B h = 34 + 12 x 10 6 M o , consistent with the value of 
35 + 8 x 1O 6 M derived via spectroscopic reverberation map- 
ping by Peterson et al. (2004). 

Already in the 1970s and 1980s numerous studies revealed 
that the Hp line profile of Arkl20 is strongly variable and ex- 
hibits two peaks in the rms spectra and a prominent red wing, 
which is neither present in Ha nor in Hy (Fig. |7). As discussed 
by Doroshenko et al. (1999), the BLR of Arkl20 may show ac- 
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Fig. 9. Cross correlations of Ark 120 light curves for B and NB, 
and for B and the synthetic UB (= NB - 0.5 V). The dotted lines 
indicate the error range (±l<x) around the cross correlation. The 
red shaded area marks the range used to calculate the lag T cen t 
(vertical dashed line). 
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Fig. 10. FR/RSS cross correlation results of Ark 120. left: B 
and NB, right: B and synthetic HB. Each histogram shows the 
distribution of T cent obtained by cross correlating 2000 flux ran- 
domised and randomly selected subset light curves. The median 
of this distribution is taken as lag T cent . The red shaded area 
marks the 68% confidence range used to calculate the errors of 
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cretion inflows in addition to the virialised motion, and therefore 
any black hole mass estimates should be considered with some 
reservation. 
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Fig. 11. B versus V fluxes of Ark 120, measured in a 7.5" aper- 
ture. The fluxes for each night (fat dots) follow a linear gradient 
yielding the range for the AGN slope. As host slope we adopted 
the range found by Sakata et al. (2010). The intersection of the 
two slopes (red shaded area) defines the host flux. The vertical 
dotted line marks the host flux found in an aperture of 5" x 7'.'6 
via modelling of HST F550M images by Bentz et al. (2009). All 
data are corrected for galactic foreground extinction. 



3.2.4. Host-subtracted AGN luminosity 

To separate AGN and host luminosity contributions, we applied 
the flux variation gradient method (Choloniewski 1981). Fig.fTTI 
shows the B versus V fluxes of Ark 120, corrected for galactic 
foreground extinction (Schlegel et al. 1998). Table [2] summa- 
rizes the results. A linear least-squares fit to the flux variations 
yields T BV = 1.19 ± 0.06, consistent with r BV = 1.20 + 0.10 
found by Winkler et al. (1992, using 20" aperture and adopting 
Eb-v = 0.1) and with Tbv = 1 22 + 0.02 found by Doroshenko 
et al. (2008, after correcting the values in their Table 3 for ex- 
tinction; note that due to a typo the value V = W.Tlmag in 
this Table should be V = \A.12mag, in order to be consistent 
with the plotted values in their Fig. 4). The range of host slopes 
is taken from Sakata et al. (2010). The intersection area of the 
AGN and host slopes defines the host contribution. Averaging 
over the intersection area yields a mean host flux of fB = 4.26 
+ 0.59 mJy and fV = 9.02 ± 0.62 mJy. Our host fluxes (in a 7/5 
aperture) are consistent, albeit a bit high, when compared with 
the values fB * 4.5mJy and fV « 9mJy (in a 15" aperture) 
from a FVG analysis of UBVRI data (Doroshenko et al. 2008). 
Subtracting our host fluxes from the total fluxes, we obtain the 
host-subtracted B and V fluxes and then by power-law interpola- 
tion the rest-frame AGN flux / 5100 a = 9.41 ± 0.86 mJy. At the 
distance of 141.8 Mpc, this yields a host-subtracted AGN lumi- 
nosity £ 510 oa = 13.27+ 1.21 x 10 43 erg/ s. The uncertainty of the 
AGN luminosity is less than 10%. 

The AGN luminosity during our campaign as derived with 
the FVG method is about 50% higher than the mean value 
.L5100A = 8.47 ± 0.81 x 10 4i erg/s derived with host galaxy 
modelling by Bentz et al. (2009) for the data of the spectro- 
scopic monitoring campaign (Peterson et al. 1998a). A detailed 
comparison of the parameters yields that our aperture area (7"5 
in diameter) is only 16% larger than that of Peterson's cam- 
paign (5" x 7'.'6), and the V band host fluxes agree within 
the uncertainties as shown in Fig. QTJ (extinction corrected: 
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fV ~ 9m Jy from our data and F550M ss 10m Jy from the 
HST/spectroscopic data base). However, the total observed (not 
extinction corrected) fluxes F 5lQ0 ^ differ significantly, being 
13.4 X 1CT 15 ergs' 1 cm~ 2 A~ l for our campaign and lying in the 
range 7.82- 10. 37xlO _15 ergs~ 1 cra~ 2 A~ 1 during Peterson's eight 
years campaign, hence are 30% - 70% larger at our campaign. 
This leads us to conclude that the AGN luminosity is ~50% 
higher during our campaign. 



4. Discussion 

The tests on PG0003+199 and Ark 120 demonstrate the feasibil- 
ity of photometric reverberation mapping to measure the BLR 
size. The case of Arkl20 demonstrates that the continuum con- 
tribution (~50%) to the narrow-band NB can successfully be 
estimated and corrected (H/3 = NB - 0.5 V), while the uncor- 
rected cross correlation (fi/NB) shows two distinct peaks, one 
peak from the autocorrelation of the continuum and one peak 
from the emission line at lag r cent consistent with the literature 
data. In the case of Arkl20, it turns out that the host contribution 
to the V band is 50% and the AGN+host continuum contribution 
to the NB is about 50%, too. Therefore, due to this coincidence, 
the V scaling factor of 0.5 works well. The general procedure, to 
recover a pure H/3 light curve, depends on the AGN/host contrast 
CAGN/host in the filters used and the AGN/line ratio CACN/ime in 
the NB. Thus, a proper treatment of the relevant filter light curves 
would be: 

• Remove the constant host emission from the light curves: To 
this end, estimate the host contribution in each filter, for in- 
stance via FVG analysis. Subtract the host contribution from 
the light curves. Now the (renormalized) light curves mea- 
sure the AGN continuum and in the NB also the line emis- 
sion. 

• Extract the line light curve from the NB light curve: To this 
end, estimate the AGN continuum contribution to the NB, 
for instance from a spectrum also taking aperture differences 
into account. Then, using a suitable filter next to the NB, e.g. 
V for H/3, subtract the AGN continuum light curve from the 
NB light curve. 

This procedure requires good signal/noise data. Applying the 
procedure on Arkl20 for a range of CAGN/host and Cacn/Hhc pa- 
rameters, we obtained cross correlation results in agreement with 
those of the straight forward method (H/3 - NB - 0.5 V). This 
indicates that the lag determination is quite robust against pa- 
rameter choice to remove the host and AGN continuum. A more 
detailed simulation of these issues would be interesting, but is 
beyond the scope of this paper. 

We now discuss the accuracies achieved for the two objects, 
the importance of good time sampling, and the improvement po- 
tential for the Rblr - L relationship. 

The BLR sizes determined by photometric monitoring agree 
within a few percent with those from spectroscopic reverberation 
mapping. The Rblr measurement accuracy we derived from the 
FR/RSS method is 12% for the short light curve of PG0003+ 199 
and 7% for the longer monitored Ark 120, while the uncertainty 
of the black hole mass (~35%) is mainly due to the adopted 25% 
uncertainty of the velocity dispersion from single epoch spectra. 

Our Rblr measurement accuracy appears exceptionally high 
compared with typical Rblr errors (~30%) reported so far from 
spectroscopic reverberation data. Adopting that the determina- 
tion of the emission line lag via spectroscopic monitoring is 
more precise than via photometric monitoring, the difference in 




median sampling [days] 

Fig. 12. Dependence of T cen , and its errors from the median sam- 
pling. 



the error ranges may be explained by two effects. Firstly, com- 
pared to our photometric campaign, most of the available spec- 
troscopic campaigns were performed with lower sampling rate, 
but over longer time spans. Thus they catch several variability 
events and implicitly determine an average lag. Then the re- 
ported errors include not only the formal measurement errors, 
but also the possible scatter of the intrinsic lags. On the other 
hand, a low sampling rate might limit the potential to reach 
small errors. Then our high Rblr measurement accuracy might 
be favoured by the high sampling rate (median 2 days) of the 
light curves. 

To check the influence of the sampling rate on the errors, 
we performed the following test. From the entire light curves of 
Arkl20 we created three regular subsets with a median sampling 
of 4, 8, and 16 days and determined for each subset T cent and its 
errors using the FR/RSS analysis. It turns out that poorer sam- 
pling results in substantially increased errors. Even worse, for 
the poorest data set (sampled with median 16 days) T cen , could 
not be recovered. The reason is that the two cross correlation 
peaks seen in Fig. [9] were not anymore separated. Therefore, 
in order to quantify the error dependence from the sampling, 
we applied a constraint in FR/RSS analysis. We excluded those 
^cent values outside the range 30-70 days. This enables us to 
determine T cen , even in the poorly sampled data set. The result 
of this procedure is shown in Fig. Q~2] The test confirms that 
the sampling rate plays a crucial role to achieve small errors. 
More important, if the continuum contributes significantly to the 
narrow-band, then photometric reverberation mapping requires 
well sampled light curves, because otherwise the cross correla- 
tion may fail to disentangle the line echo from the autocorrelated 
continuum contribution. 

We caution that the formal errors from the FR/RSS analy- 
sis do not account for potential systematic errors caused, for in- 
stance, by gaps in the observed light curves. We illustrate this 
with the light curves of Arkl20. Despite the dense median sam- 
pling of the light curves by 2 days, the light curves of Arkl20 
have a gap of 18 days in December 2009 (Fig. [8}. While the 
continuum brightening occurred in November, the line brighten- 
ing, as far as we can see it, is delayed to January 2010, yielding 
the measured lag of T cen , ~ 4&days. However, because of the 
observational gap we do not know whether the line brightening 
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Fig. 13. Artificial B and NB light curves of Arkl20, with addi- 
tional data points (black dots) filling the observational gap of 18 
days end of December 2009. The data points are chosen so that 
they mimic a strong line echo during the gap. 



occurred already in December 2009. To test the effect, if there 
were a strong line echo in December 2009, we took the B and 
NB light curves and filled the gap with the artificial bump in the 
NB light curve as shown in Fig. [13] The bump may be some- 
what strong, but it was chosen to better illustrates the potential 
effects here. One can see already in the light curves that the lag 
between the B band increase in November and the NB increase 
in December is about 30 days. The cross correlation of these 
artificial light curves shows two major peaks, one from the au- 
tocorrelated continuum at lag zero and one from the line at lag 
31 days (Fig. [Hi. The FR/RSS analysis yields T cent = 30.2+J° 
days. The difference between this lag and T cent = 48 days derived 
from the original light curves is 18 days which is just about the 
size of the gap. From this example we conclude that in case of 
a large observational gap one should take a potential systematic 
error of T cent into account, which is about the size of the gap. In 
this example the light curve is short in the sense that it covers 
only one pronounced continnum event and one subsequent line 
echo. Therefore the missing information during the gap cannot 
be filled. But in long lasting light curves which cover several 
variability events one may expect that on average the effective 
gaps are smaller. These considerations about systematic errors 
hold for both photometric and spectroscopic reverberation data. 
Therefore, they do not question the successful applicability of 
photometric reverberation mapping. 

Reverberation mapping, in general, relies on the time res- 
olution, and therefore it is vital that the light curves are suffi- 
ciently well sampled. Reminding Nyquist's Theorem, the sam- 
pling rate SR of the feature of interest, i.e. the lag, must be SR 
= t i a g ft sam pUng > 2.2. In the observational practice, data points 
are lost due to poor weather or technical problems, hence one 
should aim at SR > 3. So far, however, many BLR sizes have 
been derived from sparsely sampled light curves (SR < 3) with 
a limited potential to reach better than 30% accuracy. It could 
be that for undersampled data any quoted uncertainties of 30% 
are too optimistic. For example, the recent spectroscopic rever- 
beration campaign with improved daily sampling on NGC 405 1 
revealed that its BLR size is actually a factor ~3 smaller than 
inferred from earlier undersampled data (Denney et al. 2009b). 
Therefore, we suggest that some more of the BLR sizes mea- 
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Fig. 14. Top: Cross correlation of the complemented B and NB 
light curves of Ark 120 shown in Fig. [13] The dotted lines indi- 
cate the error range (±lcr) around the cross correlation. Bottom: 
Distribution and 68% confidence range of r ce „,, obtained by 
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sured so far can be improved with well sampled reverberation 
data. 

Applying the flux variation gradient (FVG) method to B and 
V band data and adopting a host colour range determined for 
other local AGN, we find that for both objects PG0003+ 199 and 
Arkl20 the host-subtracted AGN luminosity L 5l00 ^ differs sig- 
nificantly from previous measurements. While for PG0003+199 
we find a stronger host contribution leading to 50% lower AGN 
luminosity, for Arkl20 our host estimate is consistent with oth- 
ers but the AGN was in a 50% brighter state during our cam- 
paign. The AGN/host luminosity contrast of both objects is rel- 
atively low, LAGN/Lh osr ~ 1 for the 7" - 8" apertures used. 
The uncertainties of L 510 oa are 9% and 12% in Arkl20 and 
PG0003+199. These uncertainties include the measurement er- 
rors, the AGN variations and the uncertainty of the host flux. For 
Lagn I Li W st ~ 1 , the errors are dominated by the uncertainty of 
the host flux. It depends on the error of the AGN slope and the 
range adopted for the host slope. If improved image quality al- 
lows us to use smaller photometric apertures containing less host 
flux, the uncertainty caused by the range of host slopes may be 
reduced. Also, applying the FVG method to more than one fil- 
ter pair, for instance from UBVRI monitoring data, allows one to 
construct several independent AGN slopes and thus to determine 
the host contribution in a consistent manner, largely independent 
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Fig. 15. Rblr versus L, from Bentz et al. (2009). Shown is 
a zoomed portion containing the two sources PG0003+199 
and Arkl20 . For PG0003+199 we used R w = /Wl-37 
(Sect. I3.1.21 >. The new measurements shift both objects clearly 
outside the range of the previous errors bars and about a factor 2 
closer to the dashed line which is the fit with a slope a - 0.519 
obtained by Bentz et al.(2009). 



of the adopted range of host slopes. For powerful AGN with 
Lagn/Li W si > 2 the uncertainty of the host contribution plays 
a minor role. 

An important result from spectroscopic reverberation map- 
ping is the relationship between the BLR size and the nuclear 
luminosity Rblr 00 L a with a predicted slope a = 0.5 (Netzer 
1990; Netzer & Marziani 2010). This relationship allows one to 
derive the virial black hole mass for high redshift AGN from 
single epoch spectra by inferring Rblr from L (e.g. Vestergaard 
2002; Netzer 2003, McLure & Dunlop 2004). While early ob- 
servations indicated a = 0.6 ± 0.1 (Kaspi et al. 2000, 2005), for 
H/3 line and 5100 A luminosity a recent analysis including host 
galaxy subtraction yields a = 0.519 ± 0.063 (Bentz et al. 2009). 
Note that the slope of the current Rblr - L relationship depends 
on the adopted cosmology, because the more luminous sources 
are at higher redshift than the low-luminosity AGN. 

The current Rblr - L relationship exhibits a large scatter 
over an order of magnitude in both Rblr and L, and many ob- 
jects have large error bars (Fig. [T5l , The new photometric rever- 
beration measurements shift the position of PG0003+199 and 
Ark 120 in this diagram outside the quoted error range of the 
previous positions. Notably the shift is larger in L than in Rblr- 
The new positions are about 50% closer to the relationship fitted 
by Bentz et al. (dashed line in Fig. [TBI ). This suggests that well 
sampled reverberation data and improved host-free AGN lumi- 
nosities may significantly reduce the scatter of the Rblr - L re- 
lationship. Certainly, a true scatter may remain, simply because 
AGN are complicated objects and the Rblr determination de- 
pends on the continuum variability pattern and the geometry of 
the BLR which produces the line echo. Future data will show 
how far the scatter can be reduced. 



5. Outlook 

While spectroscopic reverberation mapping is the only way to 
explore the details of the innermost AGN structure and the ge- 
ometry of the BLR (e.g. Kollatschny 2003a, 2003b; Kollatschny 
& Zetzl 2010), the advantage of photometric reverberation map- 



ping with suitable filters is to efficiently measure BLR sizes and 
host-subtracted luminosities for large AGN and quasar samples 
- even with small telescopes. Photometric broad- and narrow- 
band monitoring of a sample of 100 quasars (V < 18 mag) can 
be performed with 1-m telescopes in an equivalent of 3 years 
observing time. 

In the case of Arkl20 the H/3 line contributes to only 50% 
of the flux in the narrow-bandpass we used. This suggests that 
photometric reverberation mapping of emission line lags works 
even for broader bands, as long as the line contributes at least 
50% to the bandpass. Also, it is worthwhile to test the method 
even for weaker emission lines contributing less than 50% to 
the bandpass. The upcoming Large Synoptic Survey Telescope 
(LSST) is equipped with six broad band filters and will discover 
thousands of variable AGN. The Ha line is shifted into the i' 
and z' bands at z s 0.16 and z ~ 0.38, respectively and may 
be strong enough, contributing about 50% to these broad bands. 
While an approximate photometric redshift may be sufficient to 
determine the filter which contains the Ha line, it is desirable 
to take a spectrum of the AGN, in order to measure an accurate 
redshift, also for luminosity determination, and the line disper- 
sion for black hole mass estimates. In order to determine the 
Ha line lag, the neighbouring filters, which are largely free of 
line emission, may be used to interpolate and remove the con- 
tinuum variations underneath the Ha line, in analogy to what we 
did for Arkl20. Recently, an interesting statistical alternative to 
our method has been proposed by Chelouche & Daniel (201 1), 
which is specifically designed for large samples of AGN light 
curves expected to be obtained with the LSST. Consider AGN 
light curves in two bandpasses, curve X (tracing the continuum 
largely free of emission lines) and curve Y (tracing the emission 
line with underlying continuum): While our approach aimed to 
remove the continuum contribution from curve Y by subtracting 
a scaled X curve, Chelouche & Daniel propose to use the light 
curves X and Y unchanged, and to subtract the autocorrelation 
of X from the cross correlation of X and Y, in order to deter- 
mine the line lag. Numerical simulations with synthetic AGN 
light curves and the treatment of four archival PG quasar light 
curves yield lags which are (in three of four cases) consistent 
with spectroscopic results. However, the reported lag uncertain- 
ties of individual AGN are large (~50%) and only for averages of 
large ensembles of several hundred AGN the obtained lags ap- 
pear satisfying. Future optimisations, combining our approach 
with that of Chelouche & Daniel, would be intriguing. 

Finally we outline a modification of the proposals by Elvis 
& Karovska (2002) and Home et al. (2003) to determine quasar 
distances from reverberation data and thus to probe dark energy. 
The luminosity difference between the open Einstein - de Sitter 
cosmology (Qm = 0.2, £2a = 0) and the concordance cosmology 
(D. M = 0.3, Q A = 0.7) is 20%-30% at redshift 0.4 < z < 0.8 
(Riess et al. 1998). From the Rblr - L relationship we conclude 
that two quasar samples with the same Rblr distribution, one 
sample at low redshift and one sample at high redshift, should 
have the same intrinsic L distribution. Thus one may constrain 
the luminosity for different cosmologies by measuring the BLR 
size and the host-subtracted brightness of quasar samples at dif- 
ferent redshift. If for objects matching in Rblr the true scatter 
in L can be reduced to about a factor 2 (= 200%), the mean lu- 
minosity of 1000 quasars may be determined with a statistical 
accuracy of 200%/ VlOOO = 6%, enabling a 3cr-5cr detection of 
a cosmological 20%-30% luminosity difference. 
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Table 1. Source parameters. Photometry is obtained using apertures of 7/8 for PG0003+ 199 and 7/5 for Arkl20. For the monitoring 
campaigns the photometry range is indicated by the double arrow. The photometry errors of 2%-A% for individual measurements 
are determined from the scatter caused by ~10 different calibration stars from Landolt (2009). To apply the flux variation gradient 
method, the photometry is corrected for galactic foreground extinction and converted to mJy, with results listed as fB ext corr and 
fV ext corr- Values used to correct for galactic extinction are: PG0003+199: A B = 0. 153 mag, A v = 0.1 18 mag; Arkl20: A v = 0.554 
mag, Ay = 0.426 mag (Schlegel et al. 1998). For the conversion from mag to mJy, zero mag fluxes of 4266.7 Jy in B and 3836.3 Jy 
in V are used. 



Object 


Redshift 


D L 


Date 


Bobs 


Vobs 


fB ex t corr 


fVexl corr 






Mpc 




mag 


mag 


mJy 


mJy 


PG0003+199 


0.0258 


112.6 


Aug. - Nov. 2009 


14.65 ^ 14.84 




5.69 <-> 6.78 










July 25*, 2010 


14.92 ± 0.037 


14.45 ± 0.022 


5.29 ± 0.18 


7.07 ± 0.14 








June 26"', 2011 


14.78 ± 0.036 


14.33 ± 0.033 


6.01 ± 0.20 


7.93 ± 0.24 








Aug 17'*, 1979 (McAlary et al. 1983) 


13.71 ±0.029 


13.56 ± 0.021 


16.10 ±0.45 


16.17 ± 0.31 


Ark 120 


0.0327 


141.8 


Oct. 2009 - March 2010 


14.11 ^ 14.30 


13.68 <-> 13.81 


13.59 ^ 16.20 


16.89 <-> 19.13 










± 0.036 


± 0.048 







Table 2. Parameters derived from the flux variation gradient method (FVG). Tbv denotes the FVG slope obtained by a linear 
least-squares fit to the varying total fB and fV fluxes plotted in a fB versus fV diagram (Figs. l6l and fTTTi. fB tota i and fV tota i refer 
to the mean of the total fB ext corr and fV ext corr flux ranges during our monitoring campaigns. Note that fV, ola i of PG0003+199 is 
not directly measured, but inferred from fB ex , corr via the FVG slope (Fig. [6). fBi ms , and fVhost are derived from the intersection 
of the AGN and host slopes (red shaded area in Figs. ISlandll 111. The range of host slopes is 0.40 < fB/fV < 0.53, adopted from 
Sakata et al. (2010). fB ACN = fB tota , - fB hos , and fV AC N = fVtotai ~ fV hos r, with uncertainty range cr ACN = (cr^, + cr 2 ho J - 5 . 
/acw((1 + z)5100A) gives the (foreground extinction corrected) rest frame 5100A flux, interpolated from fB A cN and fV A cN for an 
AGN continuum with power law spectral shape (F v oc y"). AL^gn (5 100 A) is the monochromatic rest frame 5100A luminosity. 



Object 




fBtotai 


fBhost 


f Bagn 


fV total 


fVhost 


fVAGN 


fAGN 


((1 + z)5100A) ext corr 


AL AAGN (5100A) 






mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


mJy 


10~ 15 erg s~' cm" 2 A -1 


10 43 erg s~ l 


PG0003+199 


1.20 


6.24 


2.16 


4.08 


7.95 


4.56 


3.39 


3.52 


4.06 


3.13 




±0.04 


±0.31 


±0.29 


±0.43 


±0.31 


±0.30 


±0.40 


±0.41 


±0.47 


±0.36 


Ark 120 


1.19 


15.09 


4.26 


10.82 


18.18 


9.02 


9.14 


9.41 


10.85 


13.27 




±0.06 


±0.73 


±0.59 


±0.94 


±0.57 


±0.62 


±0.84 


±0.86 


±0.99 


±1.21 
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